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I.  IWTRODUCTION 


Today  nuch  data  exists  vlth  regard  to  the  buckling  of  plates  and 

shells.  Hie  chief  conclusions  of  both  the  theoretical  and  experimental 
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Investigations  are  outlined  In  detail  In  several  papers.  '  Unfortunately » 
the  results  of  eiqiierlments  with  circular  cylindrical  shells  are  not  at 
the  present  time  consistent  vlth  theoretical  predictions.  Hie  values  of 
buckling  loads  determined  by  tests  differ  appreciably  from  the  cooqputed 
values,  no  matter  what  theory  we  use. 

Atten^ts  to  explain  these  differences  In  behavior  have  been  made  by 

various  researchers.  Donnell^  attributed  early  failure  to  Initial 

k 

deviations  from  the  geometrically  exact  shape.  Flugge  calculated  the 
stresses  caused  In  the  shell  by  the  restraint  to  e]q;>anslon  provided  by 
the  testing  machine.  Both  of  these  workers  were  able  to  demonstrate 
that  reductions  In  falling  load  would  occur  from  these  causes  but  they 
did  not  account  for  the  large  deviations  In  observed  buckle  shape  from 
those  predicted  by  theory. 

In  part,  the  discrepancy  may  be  due  to  the  fact  that  the  buckle 
pattern  observed  Is  not  of  the  type  normally  used  in  theoretical  consider- 
atlons.  It  is  most  unusual  In  practice  for  the  buckle  pattern  to  develop 
uniformly  over  the  whole  tube.  Rather  the  characteristic  diamonds  are 
formed  over  parts  of  the  tube.  A  wide  field  of  theoretical  exploration 
would  be  opened  up  If  an  attempt  were  made  to  cover  localized  bands  of 
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diamond  buckles.  To  some  extent  this  has  already  been  attempted  by  Hi  tf. 

The  purpose  of  the  study  reported  herein  was  to  check  whether  the 
Influence  of  initial  irregularities  overshadows  all  consideration  of 
elaborate  waveforms,  as  is  suggested  by  Cox.^ 

Hie  fact  that  Initial  Ingierfectlons  have  a  pronounced  effect  on  the 
value  of  the  first  buckling  load  and  the  point  at  which  It  occurs  Is 
supported  by  a  result  published  by  the  authors  In  a  previous  pajier.^ 

In  the  study  reported  In  this  paper  an  electroformed  nickel  cylinder 
was  tested  under  compression  and  the  buckle  load  history  Is  shown  In 
Curve  1  of  Fig.  1.  This  buckling  occurred  In  the  lower  20  percent  of 
the  shell.  As  a  result  of  this  we  were  able  to  make  a  closely  fitting 


heavy  collar  which  was  glued  around  the  shell  In  such  a  manner  as  to 
reduce  the  effective  length  of  the  shell.  The  compression  test  was 
carried  out  again  and  In  this  case  failure  occurred  at  a  higher  load 
than  had  previously  been  experienced.  Subsequent  measurement  showed 
that  the  buckles  occurred  at  the  most  Irregular  positions  on  the  cylinder. 
Thus,  our  experiment  demonstrated  among  other  things  that  the  load  and 
position  at  which  buckling  occurs  Is  a  function  of  geometric  imperfection. 

The  testing  of  cylindrical  shells  is  a  subject  fraught  with  diffi¬ 
culties.  Even  if  a  family  of  shells  with  prescribed  small  local  variations 
In  form  and  property  could  be  made,  a  large  number  of  tests  would  be  re¬ 
quired.  The  test  results  achieved  on  the  family  would  be  influen***  d  by 
many  factors  other  than  those  to  be  Investigated.  For  example,  the 
machine  used,  the  Individual  setup,  the  nature  of  the  end  fittings  and  so 
on,  would  Influence  the  results.  It  Is  clear  that  variation  due  to  these 
causes  must  be  eliminated  in  any  test  series.  Hence,  some  radically  new 
approach  to  testing  of  shells  was  necessary  if  the  question  posed  was  to 
be  answered  experimentally. 

We  reconsider  the  question  of  a  family  of  specimens  in  the  following 
manner:  'n*  cylindrical  shells  manufactured  by  the  same  processes  to 

identical  specification  should  have  a  random  variation  in  irregularity 
unless  some  bias  be  introduced  by  the  inspection  process.  These  'n' 
si>eclmens  could  be  placed  together  in  a  continuous  line  and  if  wr*  were 
to  examine  them  In  this  form,  the  random  variation  would  still  occur.  If 
now  they  were  Joined  together, they  would  produce  a  longer  shell  which 
would  be  characterized  by  the  same  normal  distribution  of  irregularity 
over  Its  length  as  had  previously  been  found  in  the  individual  members 
which  had  been  bonded  together  to  form  the  single  specimen-  Thus,  the 
single  specimen  has  basically  the  character  of  a  set. 

If  this  single  specimen  can  be  tested  in  an  appropriate  manner, then 
it  should  be  possible  from  the  one  test  specimen  to  obtain  the  character¬ 
istics  of  a  family  of  cylinders.  In  previous  researches  the  authors  have 
demonstrated^  that  when  the  depth  to  which  a  buckle  is  permitted  to  develop 
is  restricted  to  the  thickness  of  the  shell, then  the  buckle  process  has  no 
deleterious  effect  on  subsequent  buckle  behavior  of  the  shell.  Moreover, 
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if  a  restriction  of  this  Kind  is  made  by  means  of  some  interior  object 
lying  close  tw  the  surface  of  the  shell, it  is  possible,  with  Increasing 
load,  to  produce  a  family  of  buckles  which  will  cover  the  entire  surface 
of  the  shell.  Figure  2  shows  a  well>developed  buckle  pattern  on  a 
cylindrical  shell  tested  in  this  manner. 

It  would  seem  reasonable  under  these  circumstances  that  for  em 
initially  stress -free  cylinder  the  number  of  buckles  formed  would  be  a 
measuz^  of  the  geometric  and  material  liig>erfectlon.  Additionally,  the 
rate  of  change  in  number  of  buckles  versus  the  applied  load  would  follow 
a  normal  distribution  and  the  average  of  the  load  at  which  buckles  first 
were  realized  and  the  load  at  which  the  shell  was  completely  filled  with 
buckles  would  correspond  closely  to  the  buckling  load  of  a  perfect 

* 

I  cylinder  of  the  nominal  thickness.  If  the  contention  of  Cox  Is  correct, 

I  this  load  should  have  the  value  0.6(Et/r)  x  2nrt  approximately. 

t 

,  We  should  anticipate  that  the  better  the  manufacturlrig  techniques 

used  for  the  shells, the  more  leptokurtlc  would  the  distribution  curve  bet 
and,  of  course,  the  poorer  the  technique, the  more  platykurtlc.  (See 
Fig.  3.)  It  Is  not  unlikely,  too,  that  with  some  methods  of  manufacture 
and  with  some  materials  a  measure  of  skewness  will  be  apparent.  Tbls 
would  certainly  be  the  case,  for  example.  If  in  the  cylinder  under  test 
the  stress-strain  curve  for  the  material  Is  not  Hookean  over  the  complete 
range  of  stresses  which  result  from  the  loading  range  required. 
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II.  NATURE  OF  THE  TEST  SPECIMENS 


The  first  cylindrical  shells  used  in  this  study  were  manufactured 

by  electroforming  (Fig.  4).  A  thin  coat  of  nickel  was  deposited  on  an 

accurate  aluminum  mandrel.  The  shells  were  separated  from  the  mandrels 

by  rapidly  cooling  the  mandrels  with  liquid  nitrogen.  The  finished 

shells  were  .CX)U  Inch,  wall  thickness;  2-906  Inches, diameter;  and  8  Inches* 

length.  The  nickel  had  a  Young's  modulus  of  approximately  2k  x  10^ 

2 

Ibs/lns  ,  as  determined  In  a  separate  test. 
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Thus,  the  classical  buckling  load  of  the  specimens  was  computed 
from  the  formula 

Ft 

P  =  0.6  —  X  2nrt  =  lU55 
cr  r 


III.  METHOD  OF  TEST 


The  specimen  was  mounted  In  the  special  test  fixture,  shown  in 
Fig.  5*  It  is  clear  from  this  drawing  that  the  fixture  differs  from  the 
normal  "encastre"  type  fixture.  There  is  a  central  core  or  mandrel.  The 
cylinder  is  in  close  proximity  to  this  —  the  gap  being  .004  inch.  The 
mandrel  does  not  carry  any  direct  load  but  restricts  the  depth  of  buckle 
which  can  form.  This  loading  fixture  was  mounted  in  a  standard  60, 000-lb 
Baldwin-Hamllton  test  machine  and  load  was  applied  through  a  spherical 
seat  loading  pad  in  the  usual  manner  (Fig.  6). 


IV.  RESULTS 

Load  was  applied  to  the  first  cylindrical  shell  using  the  lowest 
rate  of  application  that  could  be  obtained  with  the  test  machine.  The 
specimen  was  continuously  watchv.u,and,  ininedlately  it  buckled,  the  leading 
process  was  discontinued.  The  shell  buckled  in  a  normal  diamond  pattern, 
and  the  buckles  occurred  in  a  small  localized  region  The  number  of 
buckles  was  counted  and  the  value  of  load  recorded  Next, the  load  was 
increased  by  50  lbs .  Loading  was  stopped  and  the  number  of  buckles  again 
counted-  This  procedure  was  repeated  until  the  buckles  had  almost  com¬ 
pletely  covered  the  surface  of  the  shell.  The  load  was  now  removed  and 
the  specimen  was  found  to  have  returned  to  its  initial  shape.  No  sign 
of  damage  or  permanent  deformation  could  be  found.  In  Table  i  the  number 
of  buckles  and  the  load  at  which  the  observation  was  made  is  recorded. 

The  results  are  graphically  portrayed  in  Figs.  7  and  8.  These  presenta¬ 
tions  show  that  the  variation  of  number  of  buckles  with  load  closely 
approximates  to  the  Gaussian  distribution,  which  is  characteristic  of  a 
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rcmdom  variation.  It  Is  seen^  too,  very  clearly  froa  Fig.  7  that  when 
the  change  In  number  of  buckles  Is  plotted  against  the  load  the  re¬ 
sulting  curve  has  zero  slope  at  a  point  vhlch  Is  very  close  to  the 
critical  buckling  load  estimated  from  the  classic  buckling  equation. 

Ihe  extremely  encouraging  result  obtained  with  the  first  cylindrical 
shell  led  us  to  believe  that  this  method  of  Investigation  was  Indeed 
fruitful.  As  a  consequence,  a  second  shell  was  prepared  In  the  same  manner 
as  the  first.  In  this  case  the  procedure  adopted  was  as  before  but  with 
addition  of  displacement  measurements.  'Die  deflection  versus  load  history 
Is  shown  In  Fig.  9  and  tabulated  In  Table  2.  Ihe  buckle  load  behavior  Is 
as  given  In  Table  3  and  Figs.  10  and  11.  Due  to  a  slight  failure  of  the 
clamping  arrangements  at  the  upper  end  of  the  cylinder  the  test  was  dis¬ 
continued  before  the  shell  was  coiiq;)letely  filled  with  buckles.  Ihe  limited 
data  obtained,  however.  Is  fully  consistent  with  that  obtained  on  the 
previous  test.  There  Is,  as  one  ml^t  expect  however,  a  variation  In  the 
value  of  a  which  defines  the  distribution  curve.  This  we  attribute  to 
the  fact  that  although  the  method  of  manufacture  of  the  two  shells  was 
the  same  the  second  shell  was  of  somewhat  higher  quality  than  the  first. 

A  third  nickel  shell  which  had  been  prepared  In  the  early  stages 
of  the  development  of  the  electroforming  technique  which  was  not  of  the 
same  quality  as  the  other  two  was  also  available.  Despite  Its  lower 
quality  we  thought  that  It  would  be  worthwhile  to  test  the  shell, and  the 
result  of  this  work  Is  presented  In  Table  h  and  Figs.  12  and  13.  The 
shell  exhibited  the  same  behavior  as  had  been  experienced  with  the  other 
two  but  in  this  case  the  nominal  distribution  curve  was  more  platykurtlc 
than  had  previously  been  the  case.  This  demonstrated  very  forcibly  to 
us  that  the  kurtosls  of  the  curve  Is  certainly  a  measure  of  the  Irregu¬ 
larity  of  the  cylinder. 

We  must  enqphaslze  that  a  coiiq>arlson  of  the  normal  distribution 
curves  which  resulted  from  the  three  tests  shows  the  following  In^rtant 
facts.  At  the  low  end  of  the  load  scale  there  Is  disagreement  between 
the  exx>erlmental  data  and  the  distribution  curves.  This  disagreement 
arises  because  the  true  distribution  curve  should  result  from  a  very  long 
cylindrical  shell  whereas  we  tested  a  short  one  in  which  end  effects 
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certainly  have  some  significance  in  relation  to  the  buckle  load  and 
number  at  the  ends  of  the  cylinder. 

The  modes  for  the  three  distribution  curves  are  located  at  1300, 

1^50  and  1^75  lbs.  The  computed  buckling  load  for  the  cylindrical  shells 
Is  1455  lbs.  Thus,  we  see  that  the  mode  value  is  in  no  case  different 
from  the  theoretical  buckling  load  by  more  than  10  percent  A  10  percent 
variation  in  this  value  corresponds  to  a  variation  of  thickness  of  5  per¬ 
cent  or  two  10 -thousandths  of  an  inch. 


V.  CONCLaSTONS 


We  have  shown  by  restricting  the  depth  to  which  a  buckle  can  develop 
that  it  is  possible  to  cause  buckling  to  occur  not  only  at  an  isolated 
point  on  the  surface  of  the  shell  but  over  ail  tne  body  The  fact  that 
the  variation  in  the  number  of  buckles  which  occur  as  a  function  of  load 
follows  the  normal  distribution  curve  nno  that  the  mean  of  the  smallest 
and  greatest  buckling  loads  approximates  closely  the  classical  critical 

buckling  load  is  clear  evidence  that  the  random  distribution  of  imper¬ 
fections  which  occur  in  cylinders  is  a  main  contributory  cause  to  the 
substantial  scatter  which  is  observed  in  cylinder  tests.  This  point  of 
view  is  substantiated  by  the  observations  recorded  in  a  previous  paper 
when  two  distinct  buckling  loads  were  obtained  from  the  same  cylinder  by 
the  simple  expedient  of  preventing  the  first  buckle  occurring  on  the 
second  test.  The  closeness  of  the  critical  buckling  load  to  the  classic 
buckling  load  seems  substantial  support  for  the  small  displacement  theory. 
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VII.  APPENDIX 


In  this  appendix  ve  repoz*t  additional  work  which  was  performed  with 
a  view  to  consolidating  that  reported  in  the  main  section  of  the  paper. 

If  variations  in  geometry  and  mechanical  property  (which  together 
constitute  irregularity) cause, as  we  have  demonstrated  them  to  do,  statis¬ 
tical  variation  In  the  number  of  buckles  observed  In  the  compressive 
buckling  of  a  cylindrical  shell  when  the  buckle  depth  Is  restricted, 
then,  by  an  argument  analogous  to  that  already  developed,  we  can  show  that 
similar  phenomena  should  occur  In  the  buckling  of  a  spherical  shell 
or  hemispherical  cap  when  we  restrict  the  depth  of  buckle.  Similarly, 
we  can  deduce  that  a  conical  frustum  tested  under  compression  with 
buckle  depth  restriction  should  not  have  a  variation  of  buckle  number 
with  load  which  follows  the  Gaussian  distribution.  This  arises  from 
the  fact  that  the  individual  frusta  which  constitute  the  frustum  tested 
do  not  have  the  same  critical  buckling  load. 

For  the  first  test  of  the  new  series  a  hemispherical  shell  was  used. 
Unfortunately,  this  hemispherical  shell  had  been  used  previously  many 
times  to  demonstrate  buckle  patterns  for  such  shells  under  external 
loading.  Thus  it  without  question  was  affected  by  the  plastic  conditions 
which  occur  at  the  crest  of  the  buckles  in  all  shells  of  single  or  double 
curvature  which  are  repeatedly  loaded  without  se'-ere  restriction  on  the 
buckle  depth  permitted.  Thus,  for  this  specimen  the  pressure  to  cause 
collapse  was  not  in  any  way  associable  with  the  critical  buckle  pres¬ 
sure  as  computed  by  any  currently  existing  theory.  On  testing  with 
buckle  restriction  and  determining  the  variation  of  the  number  of  buckles 
as  a  function  of  the  external  pressure, it  was  demonstrated  that  the  anti¬ 
cipated  Gaussian  distribution  occurs  in  this  case.  The  results  of  this 
experiment  are  given  in  Fig.  lU. 

The  second  test  made  was  on  a  conical  frustum.  In  this  case  the 
variation  in  the  number  of  buckles  as  a  function  of  load  is  as  given  in 
Fig.  15.  It  is  seen  from  this  curve  that  the  variation  is  definitely 
not  of  the  Gaussian  type. 
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It  is  well  known  that  cylindrical  shells  under  combined  loading 
behave  differently  from  cylinders  under  pure  compression.  For  exan^le, 
cylinders  under  compression  and  internal  pressure  buckle  at  a  higher 
load  than  cylinders  in  pure  compression.  The  limiting  compressive  stresst 
however,  never  exceeds  the  critical  value  as  determined  from  the  classic 
buckling  formula 


Cylinders  under  compression  and  external  pressure  buckle  at  loads  which 
are  always  lower  than  is  achieved  in  pure  compression.  Moreover, 
machined  cylinders  frequently  are  prestressed  and  then  locked-in  stresses 
are  similar  to  those  which  result  from  external  pressure.  We  conjecture 
therefore  that  a  normal  prestressed  cylinder  would  behave  in  an  analogous 
manner  to  a  cylinder  under  combined  external  pressure  and  compression; 
i.e.,  we  should  anticipate  that  such  a  cylindrical  shell  would  buckle 
in  a  similar  manner  to  a  stress-free  shell  but  that  there  would  be  a 
variation  in  the  position  of  the  mode.  To  test  this  contention,  we 
constructed  a  thin-walled  cylindrical  shell  by  machining  from  a  thick 
tube.  This  thin-walled  shell,  when  removed  from  the  mandrel  on  which 
it  was  fabricated  was  very  noticeably  oval  in  cross  section.  It  was 
forced  to  the  circular  shape  and  tested  in  the  same  way  as  before;  in 
this  case  we  obtained  a  statistical  variation  of  the  number  of  buckles 
as  a  function  of  load, but  the  mode  value  did  not  agree  with  the  critical 
buckling  load  as  determined  from  the  classic  formula  (Fig.  l6). 

We  conclude  from  the  auxiliary  tests  reported  in  this  appendix 
that  the  technique  reported  in  the  main  paper  is  sound  and  the  deduct; one 
therefrom  are  valid. 
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TAB1£  1 


Number  of  Buckles  versus  Load  for  Cylinder  No.  1 


Load  (lbs)  Number  of  Buckles 


995 

3 

1050 

11 

1100 

21 

1150 

UO 

1200 

66 

1250 

95 

1300 

161 

1350 

22U 

lUOO 

261 

IU50 

287 

10 


TABIZ  2 

Load  versus  Deflection  for  Cylinder  No.  2 


Load  (lbs) 
0 
50 
100 
150 
200 
250 
300 
350 
UOO 

U50 

500 

550 

600 

650 

700 

750 

800 

850 

900 

950 

1000 

1050 

1100 

1150 

1200 

1225 

1250 

1300 

1350 

lUOO 

1^50 


-ii 

Deflection  (10  Inches) 
0 
3 
7 
11 

15.3 

20.5 
25 

28.9 

32.2 
36 

39.8 

43. 2 

47 

50.5 

54.4 

58.2 
62  .2 
66 

70.2 
74.1 

78.8 

83 

88 

93.6 

100.1 

103.5 
117 

129.6 

146.6 
170.8 

216.2 
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TABI£  3 

Number  of  Buckles  versus  Load  for  Cylinder  No.  2 


Load  (lbs. 


Number  of  Buckles 


1225 

2 

1250 

12 

1300 

20 

1350 

38 

lUOO 

10k 

1500 

169 

TABI£  U 

Number  of  Buckles  versus  Load  for  Cylinder  No.  3 


Load  (lbs)  Number  of  Buckles 


1183 

13 

1250 

25 

1300 

37 

1350 

51 

1400 

92 

1450 

134 

1500 

168 

1550 

200 

LOAD  (lb») 


Fig.  1.  Repeated  Buckling  Behavior  of  a  Cvllndrie*! 


Fig.  2.  Phiotograph  of  Cylind^'r  with  Com|  y  Developed  Buckle  Pattern. 


15 


91 


il 


10 


total  number  of  buckles 


19 


20 


Change  in  nui 


LOAD-lbt 


Fig.  10.  Comparison  of  the  Distribution  of  the  Number  of  Buckles  with 
the  Normal  Distribution  for  the  Second  Cylindrical  Shell. 


LOAO-lbt 

Fig.  11.  Plot  of  Total  Number  of  Buckles  versus  Load  for  the  Second  Shell. 
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Fl^.  13.  Plot  of  Total  Ruaber  of  Buckles  versus 
Load  for  the  Third  Shell. 


TOTAL  NUMBER  Of  BUCKLES 


PRESSURE  (inchtt  HtO) 

Fig.  l4.  Plot  of  Cumulative  and  Rate  of  Change  Dlatrlbutlona  of  Number 
of  Buckles  versus  Load  for  Spherical  Shell. 
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RATE  Of  CHANGE  Of  NUMBER  OF  BUCKLES  PER  INCH  H^O  PRESSURE 


TOTAL  NUMBER  OF  BUCKLES 


Fig.  15.  Plot  Of  CuBulatlve  and  Rate  of  Change  Dlatrlbutlona  of  Ruaber 
of  Bucklea  versua  Load  for  Tvo  Conical  Sheila. 
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RATE  OF  CHANGE  OF  NUMBER  OF  BUCKLES  PER  POUND 


M 


LOAD  (Ibt) 


Fig.  l6.  Plot  of  Cumulative  and  Rate  of  Change  Dlatrlbutlona  of  Number 
of  Bucklea  veraua  Load  for  Streaaed  Cylindrical  Shell. 
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strong  indication  that  buckling  takes  place  first  at  the  strong  indication  that  buckling  takes  place  first  at  the 

weakest  section  of  the  shell  and  that  the  amount  of  buck-  weakest  section  of  the  shell  and  that  the  amount  of  buck¬ 
ling  is  a  function  of  the  area  of  weakness.  ling  is  a  function  of  the  area  of  weakness. 
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